Dynamics of laser-induced lattice formation in Bose-Einstein condensates 
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We present variational and numerical solutions of the generalized Gross-Pitaevskii equation for 
an atomic Bose-Einstein condensate with time-dependent, laser-induced dipole-dipole interactions 
(LIDDI). The formation of a spatially modulated ground state due to LIDDI in an elongated (cigar- 
shaped) BEC is investigated. We find that modulated densities can emerge if the laser intensity 
exceeds a critical value, which depends on the linear density of atoms in the BEC. 
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I. INTRODUCTION 



The confinement of atomic Bose-Einstein condensates 
(BECs) in optical traps [1,2] and lattices [3] is commonly 
provided by a spatial inhomogeneity of the laser intensity, 
which gives rise to gradient (or dipole) forces acting on 
single atoms. Completely different dynamical effects are 
generated in BECs by interatomic laser-induced dipole- 
dipole interactions (LIDDI) [4,5]. These effects take place 
even if the atoms are illuminated by a plane electromag- 
netic wave. One such effect, predicted by means of a 
ID variational ansatz, is that LIDDI caused by an off- 
resonant circularly polarized plane wave traveling along 
a cigar-shaped BEC results in a spatially periodic density 
modulation [5] . In the present paper we discuss how such 
density modulations can emerge from a uniform BEC 
when the laser intensity reaches sufficiently large values. 



A. Supersolids versus BECs in optical lattices 

It is tempting to refer to such a density-modulated 
BEC as a supersolid. Indeed, Ref. [5] used this term 
to emphasize the presence of both spatial modulation 
and phase coherence. Moreover, there exists also close 
connection to a BEC in an optical lattice [3]. For this 
reason we dedicate this subsection to a comparison of 
our density-modulated BEC to supersolids and BECs in 
optical lattice. 

Superfluid systems with spatially modulated densities 
have a long history. The original question, however, was 
not whether a homogeneous BEC can be spatially mod- 
ulated. Instead, the discussion started from the question 
whether a solid composed of bosonic atoms can be super- 
fluid. First, Penrose and Onsager [6] demonstrated that 
Bose-Einstein condensation is impossible in an ideal crys- 
tal. However, later it became clear [7,8] that if there were 
vacancies present in a crystal then Bose-Einstein conden- 
sation would be possible. The work of Chester [9] and 



Leggett [10] led to a deeper understanding of of Bose- 
condensed solids. However, the authors of Refs. [6-10] 
concentrated on crystals where the occupation number 
of a lattice site never exceeds 1. Thus the estimation for 
the ratio between the superfluid density and the normal 
density in such a system was rather pessimistic (~ 10~ 4 ) 
[10]. In contrast, the situation for liquid 4 He is more 
favorable. For example, Pitaevskii [11] predicted the for- 
mation of a standing wave in helium (below the critical 
temperature) flowing at a supercritical velocity in a cap- 
illary tube. The systems studied in the references men- 
tioned above are quite different from a Bose-condensed 
dilute gas. 

However, the studies of Pomeau and Rica [12] and 
Goral et. al [13] are closer to our work. Pomeau and 
Rica analyzed what could happen if the roton minimum 
in the liquid helium spectrum reaches the horizontal axis 
and even if the homogeneous system would become un- 
stable at the corresponding wavelength. Their finding is 
that as soon as the uniform ground state becomes unsta- 
ble, another ground state, spatially modulated, comes to 
existence, instead. Pomeau and Rica call such a density- 
modulated state "supersolid" , an expression already used 
by Cheng [8]. 

Goral et. al discussed quantum phases of dipolar 
bosons in an optical lattice. Their system consists of a 
dilute atomic gas where the atoms interact not only via 
the short-range s-wave scattering, but also - due to their 
large magnetic dipole moments - via long-range dipole- 
dipole interaction. This dipole-dipole interaction is some- 
what similar to the LIDDI, except that the magnetic 
dipole-dipole interaction is purely static, varies as the 
inverse distance cubed, and does not exhibit oscillatory 
(retardation) features, which are esssential in considera- 
tions of the LIDDI [4,5]. They found that in an optical 
lattice such a system can have various phases, including a 
supersolid phase. We conclude our brief discussion of su- 
persolids by mentioning that recently a supersolid phase 
of 4 He may have been observed [14]. 
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Another supcrfluid system that displays spatially mod- 
ulated densities is a BEC in an optical lattice. Here up 
to several hundred atoms can occupy one lattice site. 
As first demonstrated by Anderson and Kasevich [15], 
macroscopic quantum interference can be observed in 
such systems. As we show in the present paper, in our 
system many atoms contribute to a maximum in the den- 
sity. Therefore, we expect our system to have similar 
coherence properties as a BEC in an optical lattice. 

Due to its analogy to BECs in optical lattices we 
may call the effect considered in the present paper a 
spontaneous formation of an optical lattice in a laser- 
illuminated BEC. Such a lattice is self-sustained, since it 
is formed by interfering incident and coherently backward 
scattered light waves. 



interaction and the short-range s-wave scattering. Sec- 
tion III is devoted to a variational approach to the prob- 
lem. We assume a cigar-shaped BEC with a tight radial 
confinement and arrive at a one-dimensional equation. 
We use variational calculations to obtain approximate 
stationary solutions of this equation. In Sec. IV we solve 
the full Gross-Pitaevskii equation numerically and find 
time-dependent as well as stationary solutions. Based on 
these results we discuss a possible experiment in Sec. V. 
Furthermore, we comment on the depletion of the con- 
densate in such an experiment and mention other strate- 
gies to obtain spatially modulated densities. We also 
compare density-modulated BECs formed by LIDDI to 
coherent atom recoil lasing. We conclude by summariz- 
ing our results in Sec. VI. 



B. Formulation of the problem 



II. GROSS-PITAEVSKII EQUATION 



While Ref. [5] has suggested the intriguing possibil- 
ity of "spontaneous laser-induced periodicity" in a BEC, 
it has left open basic questions related to their observ- 
ability: (a) Is the suggested ID variational steady-state 
solution suitable for a cigar-shaped BEC, or should we ex- 
pect radial density modulation as well? (b) Is a spatially 
modulated ground state - the signature of spontaneous 
laser-induced periodicity - obtainable from a uniform ini- 
tial ground state by gradually increasing the laser inten- 
sity, even though it must compete with evaporation via 
Rayleigh scattering and possible non-adiabatic effects? 
(c) Is the this density modulation distinguishable from 
that obtained by the superradiant response of a BEC 
[16]? 

In this work we aim at a realistic investigation of laser- 
induced formation of spatially modulated densities in a 
BEC that would provide clear answers to the questions 
raised above. This is accomplished through an exhaustive 
numerical investigation that fully accounts, for the first 
time, for the 3D geometry of the cigar-shaped sample and 
the time-dependent dynamics of the formation process. 
Our numerical results show that in the proposed setting 
spontaneous laser-induced periodicity is experimentally 
feasible. They are supplemented by a simple variational 
approximation, which is shown to explain the essential 
features of the numerical results and their stability. We 
also compare and contrast the differences of the present 
LIDDI effects to superradiant light scattering by a BEC 
[16]. 



C. Outline of the paper 

Our paper is organized as follows: In Sec. II we in- 
troduce a Gross-Pitaevskii equation which takes into ac- 
count the laser-induced dipole-dipole interaction and dis- 
cuss the interplay between the long-range dipole-dipole 



The time evolution of the BEC order parameter "J" is 
governed by the 3D Gross-Pitaevskii equation 
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denotes the external trapping potential and m and a s are 
the atomic mass and the s-wave scattering length due 
to short-range interatomic forces, respectively. The last 
term in Eq. (1), which comes in addition to the terms of 
the usual Gross-Pitaevskii equation [17,18], accounts for 
the long-range LIDDI. Note that \&(r, t) is normalized to 
the number of atoms. 

For a circularly-polarized laser wave propagating along 
the z axis, the kernel of this effective interaction reads [5] 



V dd (r) = U 



(3cos 2 tf- 1) 



cos qr + qr sin qr 



(cos 2 $ + 1) 



cos qr 
qr 



(qr) 3 

cos(qr cos$). (3) 



Here, the direction of the interatomic-separation vector r 
relative to the z axis is determined by the polar and az- 
imuthal angles $ and (p, respectively, and q denotes the 
laser wave number. The characteristic energy Uo cor- 
responding to the interaction of two atoms separated by 
qr ~ 1 is proportional to the laser intensity. On the other 
hand, the rate Tn ay of the incoherent (Rayleigh) photon 
scattering per atom is also proportional to the laser inten- 
sity which provides us with the connection U = \hTn ay . 
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From the outset, we must address two questions: (i) 
Can we use in Eq. (1) the bare potential Vdd(r) of Eq. (1)? 
(ii) To what extent does the interplay between the long- 
range potential Vdd(r) and the short-range interatomic 
potential change the value of a s in Eq.(l)? 

To answer the first question we recall that according 
to Belyaev's theory of effective interactions [19], the bare 
potential should be replaced by an effective potential, 
such that the scattering amplitude produced by the effec- 
tive potential in the first Born approximation coincides 
with the exact scattering amplitude generated by the bare 
potential. We have performed such a calculation for the 
isotropic part V,(r) of the L1DD1 potential, which has the 
asymptotic forms Vi(r) w —Uo(qr)~ 2 sin(2gr) for qr 3> 1, 
and Vi(r) w - (22/15) Uo(qr)- 1 for qr < 1 [20]. Wc have 
found, that, for this isotropic part, the Born scattering 
amplitude is always close to the exact one. We conjecture 
that the same is true for the entire (anisotropic) potential 
V dd (r). 

To address the second question we first recall [21] that 
the part of Vdd being proportional to (3 cos 2 -d— l)r~ 3 has 
no significant influence on a s as long as we are far from 
the shape resonance. In contrast, Ref. [22] shows that 
the addition of the Coulomb interaction to a short-range 
potential, may dramatically change the s-wave scattering 
properties. Since for short distances the isotropic part of 
Vdd has Coulomb character, this contribution may have 
some impact on the s-wave scattering length. 

To clarify this issue, we have numerically calculated 
the phase of the scattered wave, combining the isotropic 
part of Vdd and the short-range potential. The latter 
has been taken into account by the boundary condition 
x( r Y / x( r )\r^o = — l/ a s, for the solution x(r) of the ra- 
dial Schrodingcr equation describing the scattering of two 
atoms. We have not found any dramatic change of a s 
throughout the relevant range of relative momenta of the 
colliding atoms. This result is due to the fact that the 
gravity-like part of Vdd(r) has a cutoff at r <~ which 
is much shorter than the corresponding Bohr's radius [4] 
for experimentally feasible laser intensities. Thus, wc 
may use, with reasonable accuracy, the same value of a s 
as without the laser field. 



III. VARIATIONAL CALCULATIONS 

Before we solve Eq. (1) numerically we first want to 
gain a qualitative understanding of the dynamics of the 
system described by Eq. (1). For this purpose we have 
developed a quasi- ID variational approximation which 
we discuss in this section. This approximation already 
predicts spontaneous laser-induced periodicity, that is, a 
spatially modulated ground state, when the laser inten- 
sity exceeds a threshold value. 



A. Ansatz 



We start from the fact that Eq. (1) can be derived from 
the Lagrangian L — J d 3 r£ with 
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(4) 



We assume a cigar-shaped condensate where the radial 
confinement is so tight that, regardless of the atomic den- 
sity, only the ground state of the transverse motion is 
occupied. This assumption allows us to approximate [23] 
the transverse distribution of the order parameter by 

*(r, t) « (V^r)- 1 cxphOr 2 + y 2 )/(2o- 2 )W(z, t) (5) 



with a = y/h/(mu)±). This approximation reduces the 
full Gross-Pitaevskii equation to a one-dimensional equa- 
tion for the longitudinal part i[)(z,i). 
Since for 



\z\ > q~ 



the dipole-dipole potential takes the form 
V dd {v) ^ -2U {qz)- x cos 2 (qz) 



(6) 



(7) 



we expect a spatial structure with the period ir/q to ap- 
pear in the ground state. Furthermore, we assume for 
our variational approximation that the trapping poten- 
tial Vtrap (r) does not confine the BEC in z-direction, that 
is we assume W|| w 0. We therefore pursue the ID varia- 
tional ansatz 



ip(z, t) — a(t) cos (2qz) + c(t), 



(8) 



where a(t) and c(t) are complex variables. 

To stabilize the BEC against long-wavelength pertur- 
bations [5] , we must also require that the short-range in- 
teratomic repulsion, characterized by a positive scatter- 
ing length, a s > 0, is stronger than the net laser-induced 
attraction of two atoms separated by r ^ q^ 1 - Moreover, 
we argue that this short-range repulsion term does not 
affect the dynamics at the wavelength <~ n/q, since, for 
moderate BEC densities, ir/q is smaller than the BEC 
healing length, which is determined by the BEC den- 
sity and a s . Only under these conditions the variational 
ansatz Eq. (8) describes the dynamics of the condensate 
appropriately. 
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B. Equations of motion 



C. Stationary solutions 



We now derive equations of motion for the coefficients 
a(t) and c(i). For this purpose we substitute Eq. (5) with 
Eq. (8) into the Lagrangian density Eq.(4) and perform 
the integration. Here we consider a "cigar" of size i in the 
z-direction and neglect interaction terms that do not con- 
tain the large logarithmic factor A « 2 log[^/ (2a)], which 
diverges as I — > oo. The effective Lagrangian density 



2L •* \ i ^ / * • • * \ <~» * 

— = i[c c — c c) + -(a a — a a) — \l q a a 



+g 



(c*c + \a*a) 2 + \{c*a + a* c) 2 



(9) 



contains the recoil frequency Q q = 2hq 2 /m associated 
with the momentum transfer 2%, and g = (hq)~ 1 UnA is 
the effective dipole-dipole coupling constant. The Eulcr- 
Lagrange equations corresponding to Eq. (9) read 



a 



-iQ. q a + i (g/2) (3\c\ 2 a + \a\ 2 a + cV), 



c = i (g/A) (4|c| 2 c + 3\a\ 2 c + a 2 c*) 
and thus the effective Hamiltonian density 



(10) 
(11) 



2H O I |2 
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\a\ 4 + 6 \ca\ 2 + 4|c| 4 + (ac*) 2 + (a*c)' 

(12) 



is a constant of the motion. Another dynamical invariant 
of Eqs. (fO) and (11) is the linear density of atoms 



TV : 



+ \a\ 2 /2. 



(13) 



Hence, we are left with three independent real variables, 
which can be chosen as 



n= \a\ 2 -2|c| 2 , 
s = 2V2 Re (o*c), 
p= -2v^Im(o*c). 

We then arrive at the set of equations 

n = (g/2) sp, 

S = ttgp, 

p = -s (fl q + -gn 



(14) 

(15) 
(16) 



(17) 
(18) 
(19) 



which is similar to those governing the evolution of the 
Bloch vector [24]. In particular, s is the analog of the 
Rabi frequency. A general solution of these equations 
can be found in terms of elliptic functions. 



We now focus on stationary solutions, that is fix points 
of Eqs. (17)-(19) and discuss their stability. We recognize 
that the choice s = p = is the simplest stationary 
solution. The definitions of s and p imply that either a 
or c are zero. In case a = and c^Owe find from Eq. (8) 
a spatially uniform solution. From the definition of n and 
the dynamical invariant Af we arrive at n u = —27V. The 
case c = and a^O corresponds to a purely sinusoidal 
ip(z) leading to n s = 27V. These two fix points exist for 
all values of the parameters. 

The third fix point results from the choice p = and 
s 7^ 0. This condition implies Q q +gn/2 = 0, which leads 
to n ss — —2fl q /g. Moreover, p — and s ^ demands 
a ^ and c ^ and therefore corresponds to a spatial 
modulation on top of a constant background. From the 
definitions of s and p we find 



= 2^j N 2 - (n q /gf 



(20) 



Since s has to be real, this fix point exists provided 
<?7V > fl q , that is the product of g determined by the 
laser intensity and the atom density TV exceeds a thresh- 
old value, determined by the recoil frequency Q q . 

We proceed with a standard stability analysis and find 
the perturbation eigenfrequencies 



i] u = ^Jn q (n q - gAf) , 

Vs = ^n q (fl q + gAr/2) , 



la 



( 9 N) 2 - n\ . 



(21) 
(22) 
(23) 



From the expressions for rj ss and r\ u we recognize that 
the spatially modulated solution is stable throughout its 
range of existence, whereas the uniform solution is unsta- 
ble in the same region. The uniform solution is stable if 
gAf < fl q , that is when the modulated solution does not 
exist. The sinusoidal solution is always stable. However, 
its energy determined by the Hamiltonian (12) is larger 
than the one corresponding to the other stable solution. 
Hence, the modulated solution represents the spatially 
modulated ground state of the quasi-lD BEC illuminated 
by laser light. If the condition gAf > £l q is not met, 
then the ground state corresponds to the uniform solu- 
tion. The latter becomes unstable with the increase of 
the laser intensity and/or the BEC density beyond the 
bifurcation point gAf — Q q . 



IV. NUMERICAL RESULTS 

We now turn to the numerical solution of Eq. (1), us- 
ing the second-order split-operator method. The integral 
convolution of the kernel Vdd with the matter density 
1*1 2 in the last term of Eq. (1) is calculated using the 
Fast Fourier Transform (FFT) technique [25]. 
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To circumvent the difficulties with the singularity of 
Vdd at r = we introduce a small cutoff radius r c ~ 
10 _3 g _1 . We calculate the time evolution of the con- 
densate starting from the ground state of the BEC in an 
elongated (u)j_ u>u) trap and switch on the laser inten- 
sity. In order to maximize the overlap of the initial (cigar- 
shaped) and final (spatially modulated) ground state, we 
try to keep the size of the condensate fixed when we in- 
crease the intensity of the laser field. In our numerical 
calculations we change the trapping potential by decreas- 
ing uu, so as to compensate for the laser-induced com- 
pression (electrostriction) of the condensate, see upper 
left plot in Fig. 1. Another way is to ramp-up the scat- 
tering length by approaching a Feshbach resonance at a 
rate similar to the laser switch-on. 
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thy that no oscillatory structure is numerically observed 
in the radial direction. 

We find qualitative agreement between the variational 
approach and the numerical results. Our variational 
approach predicts spontaneous laser-induced periodic- 
ity for gAf/^l q > 1, whereas from our numerical data 
we find a transition to a spatially modulated density at 
gN /fl q ~ 2.5. The agreement is not too bad consider- 
ing the fact that there are only few density peaks in the 
numerically obtained density distribution. 
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FIG. 2. Stationary solutions of Eq. (1) for a scaled laser 
intensity of 7 = 234.55 (a) and 7 = 268.06 (b). The other 
parameters are the final values of Fig. 1. Only if the laser 
intensity is large enough the density is spatially modulated. 
Again, we have used different scalings for the radial and the 
longitudinal direction. 



FIG. 1. Spontaneous laser-induced periodicity due to 
LIDDI. A BEC is illuminated by a circularly polarized laser 
beam oriented along the z-axis. The upper left plot shows 
the time dependence of the laser intensity and the longitudi- 
nal trap frequency. The time evolution of the density is shown 
in a) - e) for various times marked in the upper left plot. We 
have used scaled quantities where time is measured in periods 
of the transverse motion and lengths are measured in units 
of half a laser wavelength. The dimensionless dipole-dipole 
interaction strength 7 is defined by 7 = 2tt 2 NU / (ffi,) . The 
other parameters are TY 2 uj±/Q, q = 5 and Nqa B = 25. Here N 
is the number of particles and Q q denotes the recoil frequency 
introduced in Eq. (9). Note that we have used different scal- 
ings for the radial and the longitudinal direction. 

Figure 1 clearly shows spontaneous laser-induced peri- 
odicity: After we have increased the laser intensity and 
decreased the longitudinal frequency sufficiently enough, 
a spatially modulated density emerges. 

Moreover, we perform also imaginary time calculations 
to find the quasi-stationary density distribution for two 
different values of the laser intensity. Our numerical re- 
sults show that these density distributions are spatially 
modulated only if the laser intensity exceeds a certain 
critical threshold, see Fig. 2. For the larger laser inten- 
sity we find a quasi-stationary distribution which has a 
similar shape as the final state in Fig. le. It is notewor- 



For numerical reasons we had to limit the length of 
the cigar-shapedBEC to few laser wavelengths. We have 
not found any spatially modulated density, neither in the 
longitudinal nor in the transverse direction, if the con- 
densate is significantly wider than approximately half a 
wavelength. Our numerical results therefore do not allow 
us to answer the question whether we have to confine the 
radial direction to approximately half a laser wavelength 
in order to observe spontaneous laser-induced periodic- 
ity, or whether a sufficiently large aspect ratio can already 
lead to spontaneous laser-induced periodicity. 

V. DISCUSSION 

A possible experimental realization could involve BECs 
with a linear density of approximately 630 atoms//im. 
For 87 Rb the laser detuning may be chosen to be 2ir x 6.5 
GHz, the resonant wavelength being 0.795 /jm. A pro- 
nounced spatially modulated density should then appear 
at the laser intensity I c — 38 mW/cm 2 , that corresponds 
to Tjiay — 133 s -1 . Our variational approach [Eq.(20)] 
gives I C A = 34 mW/cm 2 , where A « 2 log[^/(2cr)]. 

We may assume exponential growth of the laser inten- 
sity at the rate of ~ 10 3 s^ 1 and simultaneous decrease 
of the longitudinal frequency luu . Under these conditions, 
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the collective oscillations in the BEC, which are excited 
due to non-adiabaticity of switching of the intensity, are 
found to be small. The total depletion of the BEC due to 
incoherent photon scattering is around 50 percent. Thus 
one can see that spontaneous laser-induced periodicity is 
possible, although accompanied by considerable losses of 
the number of condensed atoms. 



A. Quantum depletion 

We are aware of the problem of the quantum deple- 
tion of the condensate. Indeed, in a recent paper [26] 
the quantum depletion has been calculated in the Popov 
approximation for a case of a laser-radiated BEC in a 
geometry alternative to that of our present work. In con- 
trast to the present geometry, where a circularly polar- 
ized plane wave is traveling along a cigar-shaped BEC, 
Rcf. [26] considers a situation in which the wave vec- 
tor and the polarization vector of a linearly polarized 
laser beam are perpendicular and parallel to the major 
axis of a cigar-shaped BEC, respectively. As the laser- 
induced roton-likc minimum in the BEC elementary ex- 
citation spectrum tends to touch the horizontal axis, the 
quantum depletion growth rapidly, the peak being cen- 
tered at the roton-minimum position in the momentum 
space. However, in the present paper we consider, unlike 
Ref. [26], not a stationary but a time-dependent regime. 
The time of formation of the distribution of the above- 
condensate atoms in a system where both the short-range 
(s-wave scattering) and long-range (LIDDI) interactions 
are present is an important question on its own. How- 
ever, we assume that the increase of the laser intensity 
in our case is rapid enough to prevent the total depletion 
of the condensate before the spontaneous laser-induced 
periodicity, that is, a transition from a uniform BEC to 
a spatially modulated BEC, can be observed. 

Nevertheless, a certain depletion is unavoidable. Thus, 
as soon as we exceed the intensity threshold, the above- 
condensate atoms act as a seed perturbation, which will 
grow exponentially at the wave vectors where the spa- 
tially uniform solution becomes unstable. Therefore, the 
instability growth is concurrent to the emergence of the 
coherent spatial structure. Although the magnitude of 
the seed perturbation has to be calculated beyond the 
mean-field (Gross-Pitaevskii equation) approach, the in- 
stability increment can be derived from the linearized 
Gross-Pitaevskii equation [18]. In numerical simulations 
of unstable BECs the role of seed perturbation is played 
by numerical noise [27]. Stable numerical algorithms 
avoid amplification of numerical rounding errors provided 
that the numerical solution converges to the exact one, 
which is stable. However, if one reaches an instability 
range, the numerical errors begin to grow due to the 
properties of the equation being solved, and not due an 
unstable numerical algorithm. Since our numerical solu- 
tions (see Fig. 1) display a stable spatial structure, we 



may conclude that our transition from one stable branch 
(the uniform one) to another (the spatially modulated 
one) is fast enough to prevent an appreciable condensate 
depletion and the development of a dynamical instability. 

B. Alternative approach 

A better strategy to obtain similar spatially modulated 
densities is to first prepare the density-modulated ground 
state in an optical lattice formed by a far-detuned stand- 
ing wave. We then switch off the optical lattice and 
simultaneously turn on a running wave which creates 
the LIDDI. It is now this running wave that keeps the 
BEC in the initially prepared density-modulated ground 
state. Although it is technically simpler to create a spa- 
tially modulated self-sustained structure in a BEC start- 
ing from a BEC confined in an external optical lattice, 
the way studied in the present paper has its advantage. 
It clearly shows the spontaneous emergence of spatially 
modulated structures when the laser intensity is gradu- 
ally increased above a threshold value. 

C. Laser-induced lattice formation versus coherent 
atom recoil lasing 

We emphasize the key difference between the forma- 
tion of spatially modulated densities due to LIDDI in 
a laser-illuminated BEC and superradiant light scatter- 
ing by a BEC: The latter, related to coherent atom re- 
coil lasing (CARL) [16], is essentially non-stationary, re- 
sulting in the generation of phonons in a BEC. In con- 
trast, the density modulations discussed in the present 
paper are quasi-stationary. Accordingly, under the con- 
ditions of slow switch-on of the laser, we may expect a 
spatially-oscillatory ground state without triggering su- 
perradiance. On the other hand, CARL has no threshold 
in terms of the laser-beam intensity, unlike the effects 
discussed above. Hence, a weak, abruptly switched-on 
laser pulse would yield only superradiance. 

VI. CONCLUSIONS 

To conclude, our 3D numerical calculations have sub- 
stantiated the simple ID variational ansatz, predicting 
the gradual formation of an axial oscillatory ground-state 
density modulation in a laser-illuminated cigar-shaped 
BEC. Its origin can be unambiguously traced to LIDDI. 
Although the depiction of the condensate due to inco- 
herent photon scattering and due to instabilities at the 
transition from a uniform condensate to a spatially mod- 
ulated condensate is considerable, we believe that the 
creation of spatial modulations in a BEC is, as discussed 
in the present paper, experimentally feasible. We empha- 
size the similarities between the spatial modulations due 
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to LIDDI and the modulations of a BEC in an optical 
lattice. We may call the formation of a spatially modu- 
lated density a spontaneous formation of a self-sustained 
optical lattice which is formed by interfering incident and 
coherently backward scattered light waves. 
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